Hilbert-Schmidt Operators and Frames - 
Classification, Approximation by Multipliers 

and Algorithms 

Peter Balazs* 
February 2, 2008 

Abstract 

In this paper we deal with the connection of frames with the class of 
Hilbert Schmidt operators. First we give an easy criteria for operators 
being in this class using frames. It is the equivalent to the criteria 
using orthonormal bases. Then we construct Bessel sequences frames 
and Riesz bases for the class of Hilbert Schmidt operators using the 
tensor product of such sequences in the original Hilbert space. We 
investigate how the Hilbert Schmidt inner product of an arbitrary 
operator and such a rank one operator can be calculated in an efficient 
way. Finally we give an algorithm to find the best approximation, in 
the Hilbert Schmidt sense, of arbitrary matrices by operators called 
frame multipliers, which multiply the analysis coefficents by a fixed 
symbol followed by a synthesis. 
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1 Introduction 



Signal processing has become very important in today's life. Without signal 
processing methods several modern technologies would not be possible, like 
mobile phone, UMTS, xDSL or digital television. The mathematical back- 
ground for today's signal processing applications are Gahor [I9] and wavelet 
|14j theory. From practical experience it soon has become apparent that the 
concept of an orthonormal basis is not always useful. This led to the con- 
cept of frames, which was introduced by Duffin and Schaefer [16]. It was 
made popular by Daubechies [H], and today it is one of the most important 
foundations of Gabor [H] and wavelet [2] and sampling theory [T] . In signal 
processing applications frames have received more and more attention [HI EZ] ■ 

Models in physics [2J and other application areas, for example in sound 
vibration analysis [7], are mostly continuous models. A lot of problems there 
can be formulated as operator theory problems, for example in differential 
or integral equations. A interesting class of operators is the Hilbert Schmidt 
class [23] . This class forms a Hilbert space and so is similar to Euclidean space 
in the sense, that there is a way to decide if two elements are orthogonal. This 
gives some tools to work with it, as for example the orthogonal projection, 
which is a very important and easy tool for example for finding minimal 
solutions. This class of operators has therefore found also a lot of applications 
in signal processing, refer for example to [2H EB] . 

To be able to work numerically the operators have to be discretized, 
using matrices instead. Application and algorithms always work with finite 
dimensional data. So the numerical efficiency of involved algorithms is a very 
interesting questions, too. 

In this paper we will look at a connection of the theory of frames with 
the one of Hilbert Schmidt operators. It starts with prelimnaries and no- 
tations in Section [2j In Section [3] we will look at a way to describe these 
class of operators using frames. This can be done excately as in the case of 
orthonormal bases: An operator H : Ti. —>■ Ti. is Hilbert Schmidt if and only 
if ^ < oo for one (and therefore for all) frame(s) (fk)- In SectionH] 

k 

we construct Bessel sequences frames and Riesz bases for the class of Hilbert 
Schmidt operators using tensor products of such sequences in the original 
Hilbert space. We will also investigate how the Hilbert Schmidt inner prod- 
uct of an arbitrary operator and such a rank one operator can be calculated 
in an efficient way, i.e. the inner product (T, gk ® fi)'ns- This will be used in 
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Section to find the best approximation, in the Hilbert Schmidt sense, of an 
arbitrary matrix by a frame multipher given an arbitary frame. For the reg- 
ular Gabor case this has been investigated in [T7j. We will give an algorithm 
working for all frame sequences. This is an important tool for applications: 
With the Gabor multipliers an efficient way to implement time- variant filters 
have been investigated [23J. Any time- variant linear system can be mod- 
elled by a matrix, slowly time-varying linear system have a good chance 
to be close to Gabor multipliers. Other matrices might be more connected 
to a 'diagonalization' using other frames. For application irregular Gabor 
frames [22l H] , wavelet frames [H] or other analysis / synthesis systems like 
e.g. Gammatone filter banks (e.g. refer to [20]), which are mainly used for 
analysis based on the auditory system, come to mind. By the best approxi- 
mation given in this paper we offer a method and an algorithm useable for 
all systems, that form a frame, to search for the best representation of any 
matrix as a diagonal operator using these systems. We give MATLAB al- 
gorithms for this approximation. Those algorithms can be downloaded from 
http : //www . kf s . oeaw . ac . at/xxl/f initef rames/hilschmidtf ram . zip. 

2 Preliminaries and Notations 
2.1 Operators On Hilbert Spaces 

We will give only a short review, for details refer to [13]. We will denote 
infinite dimensional Hilbert spaces by Ti and their inner product with < ., . >, 
which is linear in the first coordinate. An example for a Hilbert space is the 
sequence space P consisting of all square-summable sequences in C with the 
inner product (c, d) = J^'^k ■ dk- Let f E Hi, g E Ti.2 then define the tensor 



product as an operator from 0.2 to Hi by (/ C?) g) [h) = {h, g) / for /i G 7-^2- 
2.1.1 Hilbert Schmidt operators 

A bounded operator T is called a Hilbert- Schmidt (HS) operator if there 
exists an orthonormal basis (ONB) (e„) C Hi such that 



k 
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Let HS(Hi,H2) denote the space of Hilbert Schmidt operators from Hi to 
7^2- This definition is independent of the choice of the ONB. The class of 
Hilbert-Schmidt operators is a Hilbert space of the compact operators with 
the following properties: 

• TtS(Hi,H2) is a Hilbert space with the inner product {S,T)^g = 

(Tcfc, Sek)f^, where (e^) is an arbitrary ONB. 

k 

• {f®g,h® T)^g = (/, h)^ {I, g)^ 

• \\T\\^g = ||T*||^5, and T e HS ^ T* e US. 
Furthermore it is easy [1] to show 

Corollary 2.1 Let A e HS, then {A, f g)^g = {Ag, f)^. 

For more details on this class of compact operators refer to [25] or |28j . 

2.1.2 Probenius Matrices 

Let A4rn,n bc the vector space of m x ra-matrices with the inner product, 

m— 1 n— 1 
i=0 j=0 

called the space of Frohenius matrices. Extending this definition to infinite 
dimensional matrices, we call the space of matrices, where the corresponding 
norm converges, the space of Frobenius matrices, and will denote it by 
with the norm ||.||2 2- 

Furthermore for A be a p x g and B a r x s matrix let the Kronecker 
product of A and B be the p ■ r x q ■ s matrix C with 

C'j j' I 1 I I i I ' mod r,j mod s 



A®B 



( ao,o5 
ao,iB 



ai,iB 



\ O'0,q~lB ai^q^iB 
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2.2 Frames 



The sequence {gk\k G K) is called a frame for the (separable) Hilbert space 
Ti, if constants A, B > exist, such that 



A-\\ffn<J2\(f^9k)f<B.\\f\\l yfen (1) 

k 

A is called a lower , i? a upper frame bound. The index set X will be omitted 
in the following, if no distinction is necessary. 

If a sequence (gk) fulfills only the "upper frame condition" , i.e. Yl I (/; 9k) f < 

k 

B ■ 11/11^ V / G 7^ it is called Bessel sequence. 

Let Sg : H ^ H he the frame operator Sg{f) = {fi9k) ■ Qk- S = 

k 

C*C = DD* is a positive invertible operator satisfying AIn < S < Bin and 
B^^It-c < S^^ < A'^I-j-i. If we have a frame in Ti, we can find an expansion 
of every member of 7i with this frame. Let (gk) be a frame for 7i with frame 
bounds A, B > 0. Then (gj^) = {S~^gk) is a frame with frame bounds B~^, 
A^^ > 0, the so called canonical dual frame. Every / G ?i has an expansions 
f = J2 {f,S'^9k)gk and f = J2 {f,9k)S~'^gk where both sums converge 

k(^K keK 

unconditionally in Ti. 

Let {gk} and {fk} be two sequences in H. The cross-Gram matrix Gg^j^ 
for these sequences is given by (Gg^jJ^-^ = (/m, gj), j, me K. If (gk) = (fk) 
we call this matrix the Gram matrix Gg^ . 

Using the pseudoinverse it can be shown pjj that i^j^ — ^ fk ~ 
Dj^. Furthermore it can be easily shown [1] that 

Corollary 2.2 Let (gj.) be a frame on Ti. 

^Ik ^ ^9k^9k ^ ^9k,9k^gk = ^Ik'^dk 

A complete sequence (ipk) in is called a Riesz basis if there exist con- 
stants A, B > such that the inequalities 

2 

|2 



^l|c||?< 



hold for all finite sequences (cfc). A frame ('0^) is a Riesz basis for Ti. if and 
only if there is a biorthogonal sequence (ipk), i-e. {i/jk,4'j) = ^kj for all 
This is the dual frame. 

For more details on this topic and proofs refer to [TB], [12] or [TUj . 
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3 Frames Classifiying Hilbert-Schmidt Oper- 
ators 

It can be shown very easily that 

Lemma 3.1 Let {fk) be a frame and (cj) an ONB in Ti. Let H he an operator 
H ^ H. Then 

i k i 

Proof: 

{*)>Aj2\\H*eif 
I 

and 

{*) < Bj2\\H*eif 
I 

□ 

From the proof it is clear that the right inequality is true for Bessel 
sequences, so 

Corollary 3.2 Let (fk) be a Bessel sequence and (cj) an ONB in Ti. Let H 
he an operator Ti ^Ti. Then 

E ll-^^'^llw - -^E ll-^**^*!!?!: 

k i 



As we know that an operator is Hilbert Schmidt if and only if it's adjoint 
operator is as well, we gelQ: 

Proposition 3.3 An operator H : Ti. —>■ Ti. is Hilhert Schmidt if and only if 

E w^f^wli < °° 

k 

^Thanks and credits go to Ole Christensen. In private conversation he mentioned the 
existence of the first part of this coroUary. 
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for one (and therefore for all) frame(s). 



^\\H\\ns< p^\\Hfkfn<^\\H\ 
In particular for tight frames we have 



HS 




4 Frames In The Hilbert-Schmidt Class Of 
Operators 

For two frames {fk}, {gk} an operator can be described [5] by the matrix 

This is the matrix that maps Cf^{f) i— >• Cgi^{f). It is identical to the TCS 
inner product of O and fk ® gj- This can be apphed to Hilbert-Schmidt 
operators [5] as for O G TiS the representing A4{0) is a Frobenius matrix 
with ||A<(0)||2^2 < V^^IIOII^^. This mapping : HS /(2,2) 

injective. 

Also the opposite mapping can be investigated. For a Frobenius matrix 
M define (O^*'*) (M)) h = ^ (j2 Mkj {h, tpj)^ (pk, for h e Hi then C(M) 

is in ns and \\0{M)\\^s < ^/bW \\M\\ cf. [5]. 

If the involved frames are Riesz bases these two operators are bijective. 
This can be used to investigate frames in the class of Hilbert-Schmidt oper- 
ators: 

Theorem 4.1 Let {gk)keK be a sequence in Tii , {fi)i,zj in 7^2- Then 

1. Let [gk) and (fi) be Bessel sequences with bounds B,B' , then {fi ® 
9k)ii,k)eixK is a Bessel sequence for HS {Hi, 0,2) with bound ^/ B ■ B' . 

2. Let {gk) and {fi) be frames with bounds A, B and A, B' . Then {fi ® 



9k)(i,k)&ixK is a frame forHS{Hi,H2) with bounds \/ A - A! and \J B ■ B. 
A dual frame is {fi ® g^.)- 
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3. Let (gk) and (fi) be Riesz bases. Then {fi® gk)(i^k)<^ixK is a Riesz basis 
for Ti.S(Ti.i,Ti.2) ■ The biorthogonal sequence is {fi g^). 



Proof: Suppose the operator O G HS, then 
Using the dual frames we know [5j that 



HS 



101 



ns 



< 



1 



ns y/AA' 



HS 



as the reciprocal of the lower frame bound is an upper frame bound of the 
dual frame. Therefore 



AA' \\0\\l^g < ||m(^**)(0)||^^ < BB' \\0 



\ns 



This is equal to 



AA' \\0\\l,, < 5^ 1(0, h ® gXf < BB' \\0\\l,s 



k,j 



If both sequences are Riesz bases, then 



n 



HS 

This means that [gk® fj) is a Riesz basis 



fki , fi 



fc2 ■ {9j2^9jl) 

H 



n 



J1J2- 



□ 



4.1 Efficient Calculation of the TCS Inner Product 

Let us look at the finite dimensional space C". For matrix representation 
of operators [5j in HS, the inner product {T,gk ® fi)'j-is becomes important. 
The diagonal version {T,gk ^ fk)-^s P^^Y^ essential role for frame mul- 
tipliers, see Section |5l There are several ways to calculate this HS inner 
product, which we will list in Theorem 14.51 We will first collect the following 
properties for the proof of this theorem. 

Note that with \_x\ we describe the biggest integer smaller than x. 
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Lemma 4.2 For M G M.m,n, o Frobenius matrix of dimension m x n, let 
Oec(") (M), = M^^^^^^^,^fork = 0,...,m.n~l 

With this function this space is isomorphic to C™^" with the standard inner 
product. The inverse of this function is 

Proof: The function oec^"^ is clearly linear and inverse to OTtatn- □ 
The function oec*-""^ joins the columns together to a vector. The function 
X I— OJlatn(x) separates this vector again. !Hatn(x) is well known in signal 
processing, it is called Polyphase representation there. 

We will use the standard notation for complexity of a formula 0|21]. 
writing f{x) G 0{n) if the algorithm for f{x) has order of n time complexity. 
We will abuse the standard notation a bit, as we will not only include the 
highest power term. 

Lemma 4.3 The complexity of the calculation of the following terms is 

1. inner product: Let x,y E C^, then 

{x,y) eOi2p-l) 

2. matrix-vector multiplication: Let A G J^p,q, x G C'^, then 

A-xeO{p-{2q-l)) 

3. matrix-matrix multiplication: Let A G M.p,q, B G M.q^r, then 

A- B eO{p-r ■ (2q-l)) 

4- Kronecker product of matrices: Let A G Aip^q, B G Air,s, then 

A® B eO{p- q-r ■ s) 

Proof: This can be easily be deduced from the basic definitions of these 
operators. □ 



9 



Lemma 4.4 Let A e Mr,s, B e Mp^q and C e Mq^r- Then 
{A^ ® B) ■ (occ^'^^C) = oec^^) {B ■ C ■ A) 

Proof: 

{{A- » B) (t,.c(')C)), = "f: {A^ » Bl. (.ecOC)^, 

j=0 



Let ji — j mod g and j2 - 



, so 



q-l s-1 



Jl=0j2=0 



s-1 



j2=0 



□ 



Theorem 4.5 Lei (/i;)^o a frame in C", {gk)j^o -^^^ T be a linear 

operator T : C"" —>■ C™. T/ien 

1. (T, ® 7^,)^^ = (t,ec(-) (T) , t)ec(") {gk 70>c^.„ e C» ((3mn + m - 1)) 
/or eac/i single pair {I, k). 

2. {T, gk <S> hi)^g — {Thi, gk)cm £ O {{2mn + m — 1)) for each single pair 
{l,k). 

3. {T,gk®hi)^^ = {Cg,-T-Dh^\^ G O {{L {2mn - m + 2mK - K))) 
for all values {I, k). 

I {T,gk(^hi)^s = K®^/>ec(T) G 0{{KL-{3mn-im for all 
values {I, k). 
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Proof: We will use Lemma [4.31 extensively: 

1. ) Qk hi E O (m ■ n)\ oec'-"'' is only a reordering. The complex conju- 
gation G O (m). The inner product G (9 (2 ■ (mn) — 1). So the overall sum 
is 3mn + m — 1. 

2. ) Thi E O {m (2n — 1)). The inner product G O (2m — 1). The sum is 
2mn + m — 1. 

3. ) To understand that {T,gk hi)^g = {Cg^ ■ T ■ Dj^)^ ^ use Corollary 
12.11 and the facts that 



T- 



9i 92 



9m 



T91 Tg2 



T9 



M 



and 



h2 



h 



N 



9i 92 



9m 



J 



\ 



J 



^9ni hm) 



Furthermore T- D eO {mL{2n - 1)), C ■ {TD) G O {KL{2m - 1)), so alto- 
gether we get as a sum mL2n—mL+KL2m—KL = L (2mn — m + 2'mK — K). 

4.) Using Lemma [4.41 we know that this equality is true. For the calcu- 
lation of ® C*/; we need G O [KmLn). And for the matrix vector multi- 
plication in C™" G O {KL{2mn - 1)). So overall KL{2mn - 1) + KmLn = 
KL- {3mn-l). □ 

So overall if we have to calculate the inner products for all pairs (A;, I) the 
third method is the fastest (except when n is very big and m and K very 
small). If we need only the diagonal part k = I, the second one is the most 
efficient as for using the third method we would still have to calculate the 
whole matrix and then use its trace. 



5 Best Approximation of Hilbert Schmidt Op- 
erators by Frame Multipliers 

5.1 Best Approximation By Frame Sequences 

In applications we very often have the problem to find an approximation of 
a certain object. In a Hilbert space setting, where the interesting objects 
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are in a space spanned by a frame sequence, it is known [12] that the best 
approximation is just the orthogonal projection on this space, which is in 
this case given by 

p{f) = J2{fr9k)9k = {*) 

k 

A disadvantage of this formula for practical solutions is that the dual 
frame has to be calculated. This can be time-consuming and is not needed 
per-se. But we can use the formulas we have established previously to get 

(*) = D,,C-,, = D,,Dl = D,,GlC,, 

Here we can avoid calculating the dual frame directly, instead using the 
existing algorithm to calculate the pseudoinverse. 

Theorem 5.1 Let {gk) be a frame sequence in Ti. Let V = span[gk). The 
best approximation of an arbitrary element f ETi is 

5.2 Frame Multiplier 

Frame Multipliers [6] are operators which can be described in this sense by 
diagonal matrices, i.e. operators of the form 

= ^mk {f,gk)fk, 

k 

where (m^) is a bounded sequence called the upper symbol. These are op- 
erators spanned by (7^ ® ^fc). We know from Theorem 14.11 for the Hilbert- 
Schmidt class of operators 

1. that this is a Bessel sequence for Bessel sequences (^r;.) and (fk) 

2. and because every sub-family of a Riesz basis is a Riesz sequence, that 
for Riesz bases, (7^ ® gk) is a Riesz sequence. 
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5.2.1 The Identity As Multiplier 



Let us state an interesting question regarding multipliers, which can be 
adressed by the approximation algorithm presented later. Can the identity 
be described as multiplier? It can be easily shown that 

Lemma 5.2 If and only if the identity is a multiplier for the Bessel sequence 
{gk} with constant symbol c 7^ 0, then {gk} is a tight frame. 

In the case of regular well-balanced Gabor frames, it is shown in [15] that 
if the identity can be written as multiplier for the frame {gk}-, its symbol is 
a constant sequence. In this case Lemma 15.21 is clearly equivalent to: The 
identity is a multiplier if and only if {g^} is a tight frame. We will see some 
examples later, that this is certainly not the case for general frames. 

5.3 The Lower Symbol 

For Bessel sequences {gk) and {fk) the family {g^ ® fk) is again a Bessel 
sequence. So the synthesis operator is well defined: 

C^^^j^ -.US-. I' with C{T) = (T, gk ® J,)^^ 

Using Theorem 14.51 we can express that inner product as 

C{T) = {T,gk0j,)^^ = {Tj,.,gk)^ 

So we define 

Definition 5.1 Let {gk) and (/^) be Bessel sequences for Ti, then the lower 
symbol of an operator T e HS is defined as 

(Tl{T) = {Tfk,gk)n 

The function ai '■ HS is just the synthesis operator of the Bessel 
sequence gk ® fk in T^S and therefore well defined in The name is used in 
acoordance with the use for Gabor multipliers, which are frame multipliers 
for Gabor frames, refer for example to [18]. 

Let us suppose that these elements {gk®fk) are not only a Bessel sequence, 
but even a frame sequence. This has to be supposed and cannot not be 
deduced from the status of {fk) and {gk) of being frames, as Theorem 14. II is 
of no use here. Subsequences of frames don't have to be frames again [T^ . 
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In this case the best approximation can be found by using the analysis 
and the dual synthesis operator for the projection on the closed span of the 
elements V = span | ® fk}, which are exactly those TCS operators that 
can be expressed as frame multipliers with the given frames. Let Qk be the 
canonical dual frame of ® fj^ in V then the best approximation is 

Pv{T) = J2{T,gk0j,)^^Qk = J2^LiT)Qk 

k k 

As {T,gk (S> fk)-j-is j^^^ analysis coefficients with a frame sequence, we 
know ||cT(7||2 < II'^l(2^)II2 for any other coefficients au such that the projection 
Py can be expressed in this way and hence the name 'lower symbol'. Also 
for bounded operators T which are not in TiS this inner product is defined 
and bounded by \\T\\q,^\/BB'. 

5.3.1 Approximation Of Matrices By Frame Multipliers 

In infinite-dimensional spaces not every subsequence of a frame is a frame se- 
quence, but in the finite-dimensional case, all sequences are frame sequences. 
So we can use the ideas in the Sections 15.11 and 15.31 and apply it to frame 
multipliers. All codes in here can be found online at 

http : / /www . kf s . oeaw . ac . at/xxl/f initef rames/hilschmidtf ram . zip. 

We want to find the best approximation (in the Frobenius norm) of a 
m X n matrix T by a frame multiplier with the frames {gk)f=i ^ and 
(/fc)fcLi ^ C™"- This whole section is a generalization of the ideas in [17] . 

Algorithm: 

1. Inputs: T, D, Ds 

T is a m X n matrix, D is the n x K synthesis matrix of the frame 
(gk)- This means [3] that the elements of the frame are the columns of 
D. Ds is the synthesis matrix of the frame (fk)- Often we will use the 
case {fk) = {9k) so Dg = D by default. 

2. Lower Symbol : 

Using Lemma 14.31 the most efficient way to calculate the inner prod- 
uct {T,gk fk)'Hs (^/fei fi'A;)c"- T^^^^ can be implemented effectively 
using the matrix multiplication by 

(MATLAB :) lowsym(i) = conj (D( : , i) ' * (T*Ds ( : , i) ) ) ; 
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3. Hilbert Schmidt Gram Matrix : 

We calculate the Gram matrix of {gu® fk) 

i.Gns)i^k = {9k®lk.9i®li)ns = idk, 9i)n- {fiJk)n = (^aJ/.fc- (G'/J^ 

(MATLAB :) Gram = (D'*D) .*((Ds'*Ds) . O ; 

If (gk) = ifk) then 

{Gns)i^k = \ {9k,gi)nf 

It is more efficient to use this formula in 
(MATLAB :) Gram = abs((D'*D)) . ~2; 
as this has complexity, using Lemma 14.31 ~ ■ (n^ + 2) compared to 
the original calculation with ~ ■ {n^ + m? + 1). 

4. Upper Symbol : 

Using Theorem 15. II we get the coefficients of the approximation by using 
the pseudoinverse of the Gram matrix. In the case of frame multipliers 
the coefficients are an upper symbol a. 
(MATLAB :) uppsym = pinv(Gram) *lowsym; 

5. Outputs: TA, uppsym 

For the calculation of the approximation we just have to apply the 
synthesis operator of the sequence {gk® fk) to the upper symbol. 

K 

ik = Pv{T) = Y,^kgk®lk 

k=l 

The matrix of the operator Qk fk can just be calculated [3j by {gk)i ■ 
{fk)j- 

(MATLAB :) P = D( : , i) *Ds ( : , i) ' ; 

For an implementation of this algorithm in MATLAB see the file ApprFramMult . m. 
Example 5.1 : 

We will look at two simple example in C^. 
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1. Let A = ^ Q 5 / This is clearly a multiplier for the standard or- 

thonormal basis of C^. The sequence /i = (^\, i f2 = 

is also an ONB. But the best approximation of A with this basis is 

PvW = f l llTr. )■ So this ,s an example that not even 



^ 0.4330 4.2500 ^ 
for ONBs a frame multiplier for one basis stays a frame multiplier for 
another one. 

2. Let T = /(i(p2. and let D = ( ''^'^^['^^J ] ]■ This is a frame 

^ \^ sm(30°) 1 -1 J 

with bounds A = 0.5453, B = 3.4547 and therefore not tight. Still the 

identity can be approximated perfectly (up to numerical precision) with 

the coefficients a = (3.1547,-1.3660,1.5774). So this is an example, 

where the identity is a frame multiplier for a non-tight system. 

The MATLAB-codes for these examples can be found in testapprf rammult .m. 
Example 5.2 : 

We will now use this algorithm for the approximation by regular Gabor 
multipliers as presented in [17|. We are using a Gauss window in with 
n = 32. We are changing the lattice parameters a and b. The resulting 
approximation of the identity can be found in Figure [H 

1. {g,a = 2,6 = 2). This is nearly a tight frame with the lower frame 
bound A = 7.99989 and the upper frame bound B = 8.00011. As 
expected the identity is approximated very well. 

2. {g, a = 4,b = 4) : This frame is not tight anymore, as ^ = 1.66925 and 
B = 2.36068 and we can see that the approximation is deviating from 
identity. 

3. {g,a = 8,b = 8) : This is not a frame anymore, but a Bessel sequence 
with B = 1.18034. At least some of the structure (the diagonal domi- 
nance) is still kept. 

4. {g, a = 16, 6 = 16) : This is not a frame anymore, but a Bessel sequence 
with B = 1.00001. All structure is (more or less) lost. 



16 



Figure 1: Using the algorithm for approximation with frame multipliers in the 
Gabor case: Approximating the identity by Gabor multiplier with Gauss window 
(n = 2>2) and changing lattice parameters. Top Left: (a = 2,6 = 2), Top Right: 
(a = 4,b = 4), Bottom Left: (a = 8,6 = 8), Bottom Right: (a = 16,6 = 16) 

This algorithm is not very efficient for the Gabor the special struc- 

ture is not used. For the regular case the algorithm presented in ^17j is 
preferable. The MATLAB-codes for these examples can be found in the file 
testapprGabmult.m. 

6 Perspectives 

The efficient calculation of the best approximation, in the Hilbert- Schmidt 
sense, of an arbitrary matrix by a frame multiplier is important for improving 
the numerical aspects of a frame, e.g.. quotient of the frame bounds, by using 
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weights. Refer to [S], where weighted frames are introduced. This connection 
is worth to be estabhshed more in the future. 

The idea for the approximation of Hilbert Schmidt operators can be ap- 
phed for irregular Gabor multiphers. Even if a lot of the structure is missing, 
some special properties can be exploited and the algorithms can be sped up. 
For first ideas refer to [1]. 
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